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@ Method and apparatus for simulating color print. 

@ There is provided a model where the illuminance of reflected light from the print is expressed with a linear 
combination of the illuminance of specular reflection light and that of internal reflection light. On the basis of the 
model, three methods are applicable to compute the colors of a print disposed in the three-dimensional space. 
In the first method, a specular reflection coefficient and an internal reflection coefficient, which are depending 
upon the wavelength, are interpolated by an angle of reflection d and an angle of deviation p, and the 
illuminance spectrum of the reflected light is subsequently determined (step 32),' Tristimulus values X, Y, and Z 
are then determined by integration of the illuminance spectrum according to the color matching functions (step 
S3). In the second method, the tristimulus values are determined for the specular reflection light, the internal 
reflection light, and the environmental light, respectively, and then interpolated by the reflection angle d and the 
deviation angle p (step 84). In the third method, the chromaticity coordinates are determined respectively for the 
specular reflection light, the internal reflection light, and the environmental light, and then interpolated by the 
reflection angle 6 and the deviation angle p (step S5). 
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BACKGROUND OF THE INVENTION 
Field of the Invention 

5 The present invention relates to a method of and an apparatus for simulating a color print with an output 

device, such as a display device and a printer. 

Description of the Related Art 

10 In a reproduction process of a color print with various output devices like a display device and a printer, 
it is preferable to simulate colors as close as possible to those of the actual print. There are some known 
method of color reproduction or color simulation such as Murray-Davis's equation, Yule-Nielsen's equation, 
and Neugebauer's equation. The Neugebauer's equation converts dot percents of four color separations Y, 
M, C, and K into color components R, G, and B. 

76 The above conversion equations are, however, established on the basis of idealized models and in 
many cases, can not reproduce the actual colors in a printed image. The conventional models are often not 
applicable to color simulation under the condition of observing a printed image disposed in the three- 
dimensional space. 

20 SUMMARY OF THE INVENTION 

An object of the present invention is thus to provide a method of and an apparatus for simulating colors 
of a printed color print image disposed in the three-dimensional space more precisely than the conventional 
methods. 

25 The present invention is directed to a method of simulating a color of each reflection point in a color 
print observed by an observer with an output device, the color print being disposed in a three-dimensional 
space and irradiated by a light source, the light source emitting a light beam having a luminance spectrum 
4f(\) where X denotes a wavelength of the light beam. The method comprises the steps of: determining an 
illuminance spectrum \{0,p,\) of reflected light reflected at the reflection point according to an Equation: 

30 

I(0,P,X) = {Ss(X)*COS"p + Sb(X)*COS0}</>(X) + le(X) 

wherein Sb(X) and Ss(X) denote first and second reflection coefficients, respectively, p denotes an angle 
between a first direction of the reflected light and a second direction from the reflection point to the 

35 observer, n is a constant, 9 denotes an angle of reflection of the reflected light, and le(X) denotes an 
illuminance spectrum of environmental light observed by the observer; integrating the illuminance spectrum 
\(0,p,\) over a wavelength range of visible light with respective three color matching functions of a device- 
independent colorimetric system, to thereby obtain tristimulus values in the device-independent colorimetric 
system representing the color of the reflection point; and producing an image of the color print with the 

40 output device as a function of the tristimulus values. 

In a preferred embodiment of the present invention, the first reflection coefficient Ss(X) regarding 
specular reflection light depends upon the wavelength X. 

The method further comprises the step of: converting the tristimulus values to color data in a specific 
colorimetric system of the output device. 

45 The illuminance spectrum le(X) of environmental light is equal to zero. 

According to an aspect of the invention, a method comprises the steps of: determining target tristimulus 
values X, Y, and Z of a device-independent colorimetric system for reflected light reflected at the reflection 
point according to Equations: 

50 X = k{cos"p • Xs + cose • Xb + Xe} 

Y = k{cos V • Ys + cose • Yb + Ye} 
Z = k{cos"p • Zs + cosd • Zb + Ze} 

where (Xs,Ys,Zs), (Xb.Yb,Zb), and (Xe,Ye,Ze) respectively denote first, second, and third sets of tristimulus 
55 values, k is a first constant, p denotes an angle between a first direction of the reflected light and a second 
direction from the reflection point to the observer, n is a second constant, and 6 denotes an angle of 
reflection of the reflected light; and producing an image of the color print with the output device as a 
function of the target tristimulus values X, Y, and Z. 

3 
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In a preferred embodiment, the method further comprises the step of: converting the target tristimulus 
values X. Y, and Z to color data of a specific colorimetric system of the output device. 

The third set of tristimulus values (Xe,Ye,Ze) regarding environmental light are equal to zero. 

According to another aspect of the present invention, a method comprises the steps of: determining 
5 target chromaticity coordinates (x.y) of a device-independent colorimetric system for reflected light reflected 
at the reflection point according to Equations: 

(x.y) = Cs- (xs,ys) + Cb- (Xb.yb) + Ce^ (Xe.ye) 

10 

^ cos'» p -Ts 

V^S 

COS" p Ts + cos e Tb + Te 

15 



20 



Cb = 



cos H • Tb 



cos" p -Ts + cos e Tb + Te 



— — " 

cos" p Ts + cos 8 Tb + Te 

25 

where (Xs,ys). (Xb,yb). and (Xg.ye) denote first, second, and third chromaticity coordinates, Ts, Tb. and Tg are 
constants, p denotes an angle between a first direction of the reflected light and a second direction from the 
reflection point to the observer, n is a constant, and e denotes an angle of reflection of the reflected light; 
30 and producing an image of the color print with the output device as a function of the target chromaticity 
coordinates (x,y). 

The method further comprises the step of: converting the target chromaticity coordinates (x,y) to color 
data of a specific colorimetric system of the output device. 

The third set of chromaticity coordinates (Xg.ye) and the constant Te regarding environmental light are 
35 equal to zero. 

The present invention is further directed to an apparatus for simulating a color of each reflection point in 
a color print observed by an observer, the color print being disposed in a three-dimensional space and 
irradiated by a light source, the light source emitting a light beam having a luminance spectrum ^(X) where 
X denotes a wavelength of the light beam. The apparatus comprises: means for determining an illuminance 
40 spectrum l(0.p,X) of reflected light reflected at the reflection point according to an Equation: 

l(0,p,x) = {Ss(x).cos"p + Sb(x).cosa}<^(x) + le(X) 

wherein Sb(X) and Ss(X) denote first and second reflection coefficients, respectively, p denotes an angle 
45 between a first direction of the reflected light and a second direction from the reflection point to the 
observer, n Is a constant, d denotes an angle of reflection, and le(X) denotes an illuminance spectrum of 
environmental light observed by the observer; means for integrating the illuminance spectrum 1(0, p,X) of the 
reflected light over a wavelength range of visible light with three color matching functions, respectively, of a 
device-independent colorimetric system, to thereby obtain tristimulus values in the device-Independent 
50 colorimetric system; and output means for producing an image of the color print as a function of the 
tristimulus values. 

The present invention is still further directed to an apparatus comprising: means for determining target 
tristimulus values X, Y, and Z of a device-Independent colorimetric system for reflected light reflected at the 
reflection point according to Equations: 

55 

X = k{cos"p • Xs + cosd • Xb + Xe} 
Y = k{cos"p • Ys + cose • Yb + Ye} 
Z = k{cos"p • Zs + cose • Zb + Ze} 
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where (Xs.Ys.Zs), (Xb.Yb.Zb), and (Xe.Ye.Ze) respectively denote first, second, and third sets of tristimulus 
values, k is a first constant, p denotes an angle between a first direction of the reflected light and a second 
direction from the reflection point to the observer, n is a second constant, and 9 denotes an angle of 
5 reflection of the reflected light; and output means for producing an image of the color print as a function of 
the target tristimulus values X, Y, and Z. 

The present invention is further directed to an apparatus comprising: means for determining target 
chromaticity coordinates (x,y) of a device-independent colorimetric system for reflected light according to 
Equations: 

70 

(x,y) = Cs^ (xs.ys) + Cb* (Xb.yb) + Ce» (Xe.ye) 



75 



cos" p -Ts 
COS" p Ts + cos 9 Tb + Te 



20 



25 



30 



40 



^ COS H • Tb 

Cb = 

cos" p 'Ts + cos e Tb + Tc 



cos" p -Ts + cos e Tb + Te 



where (Xs.ys), (Xb.Vb). and (Xe,ye) denote first, second, and third chromaticity coordinates, Tg, Tb, and T© are 
constants, p denotes an angle between a first direction of the reflected light and a second direction from the 
reflection point to the observer, n is a constant, and d denotes an angle of reflection of the reflected light; 
and output means for producing an Image of the color print as a function of the target chromaticity 
35 coordinates (x,y). 

These and other objects, features, aspects, and advantages of the present invention will become more 
apparent from the following detailed description of preferred embodiments with the accompanying drawings. 



BRIEF DESCRIPTION OF THE DRAWINGS 



Fig. 1 shows reflection of incident light by ink on a printing paper; 

Fig. 2 shows the relationship between a light source, a print, and a vision system (or an observer); 

Fig. 3 is a chromaticity diagram of the chromaticity point Ps(xs,ys) of the specular reflection light Ls and 

the same Pb(Xb.yb) of the internal reflection light Lb under the ideal condition with no environmental light; 
45 Fig. 4 conceptually shows a measurement system for reflection properties of a printed image; 

Fig. 5 is a graph showing an example of the illuminance spectra of the reflected light measured; 

Fig. 6 is a graph showing the chromaticity data of the reflected light obtained for coated paper samples; 

Fig. 7 is a graph showing the chromaticity data of the reflected light obtained for art paper samples; 

Fig. 8 is a block diagram showing a computer system embodying the invention, which is used for 
50 reproducing a color print disposed in the three-dimensional space; 

Fig. 9 shows a three-dimensional orientation in color reproduction of a color print used in a ray tracing 

technique; and 

Fig. 10 is a flowchart showing a procedure for computing color data (RGB data) by the ray tracing 
technique. 

55 
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DESCRIPTION OF THE PREFERRED EMBODIMENT 



The description below regards: three methods of computing colors from reflection characteristics of a 
color print; a method of measuring illuminance spectra of the reflected light, or illuminance of the reflected 
5 light with respect to each wavelength; and a method of practical color reproduction. 

A. First Method of Color Computation 

Fig. 1 shows reflection of incident light by ink on a printing paper. An incident light beam Li entering the 
70 ink is reflected via two different paths as shown in Fig. 1. A first one is specular reflection light Ls, which is 
a reflection of the incident light Li from a boundary between the surface of ink and the air, A second one is 
internal reflection light Lb. A light beam passing through the surface of ink is scattered by pigment particles 
of the ink and emitted as the internal reflection light Lb. 

Fig. 2 shows the relationship between a light source, a print, and a vision system (or an observer). The 
75 position of the vision system is generally deviated from the direction of the specular reflection light Ls by an 
angle of p (hereinafter referred to as the angle of deviation p) as clearly seen in Fig, 2. The line drawn 
between a reflection point on the surface of the print and the vision system, which corresponds to the 
direction of observation, may not exist on the plane defined by the incident light beam Li and the specular 
reflection light Ls. In such a case, the angle of deviation p is defined by the angle of the specular reflection 
20 light Ls and the direction of observation. 

The illuminance of the reflected light observed by the observer in Fig. 2 is a linear combination of the 
illuminance of the specular reflection light Ls and that of the internal reflection light Lb. which is expressed 
as: 



26 l(0.pA) = ls(0.P.X) + lb(0,P.X) (1) 

where X, \(e,f>,\), \s{e,p,\), and lb(0,p,X) respectively denote a wavelength of light, an illuminance spectrum or 
a spectral illuminance of the reflected light observed, an illuminance spectrum of the specular reflection 
light Ls, and an illuminance spectrum of the internal reflection light Lb. The illuminance spectra l(0,p,X), U- 
30 (0,p.X), and \b(e,p,\) depend upon the angle of reflection d, the angle of deviation p, and the wavelength X. 

On the assumption that the illuminance spectra ls(e,p.X) and Ib(^.p.X) are respectively divided into angle 
components and wavelength components, Equation 1 is rewritten as: 



l(d.p.X) = {Ss(X). As(0.p) + Sb(X). Ab(e.p)}<»(X) (2) 

where Sb(X) denotes a coefficient of internal reflection, Ss(X) denotes a coefficient of specular reflection, <^>(X) 
denotes a luminance spectrum or a spectral luminance depending upon the wavelength X of the incident 
light, As(d,p) denotes an angle component of the specular reflection light Ls, and Ab(0.p) denotes an angle 
component of the internal reflection light Lb. 

40 The illuminance of the specular reflection light Ls has the maximum value at the point of observation 
which is in the direction of the reflected light with the angle $, and abruptly decreases in response to the 
deviation from the direction. Referring to such phenomena, the angle component A^(e,p) of the specular 
reflection light Ls in Equation 2 may be expressed as an n-th power of cosp. where n represents a constant 
experimentally determined, while the angle component A^{e,p) is expressed as cose. Upon such assump- 

45 tions. Equation 2 is further rewritten as: 



l(e.pA) = {Ss(X).cos"p + Sb(X)«cos0}<»(X) (3) 

In actual observation of a printed image, there exist, other than a standard light source, natural light and 
50 Other environmental light, which come to the observation point directly or after being reflected from the 
printed image. A general expression of the illuminance spectrum 1(6, p,\) of the observed reflected light is 
given by addition of an illuminance spectrum le(X) of the environmental light to the right hand side of 
Equation 3 as follows: 



55 1(0. P. X) = {Ss(X).COS"p + Sb(X)*COS0}</)(X) + le(X) (4) 

In Equation 4. the coefficients of reflection Ss(X) and Sb(X) and the exponent n can be determined by 
measuring a illuminance spectrum of the reflected light with a light source which has the luminance 
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spectrum <^(X) and by analyzing the results of the nneasurement. Details of the measurement process will be 
described later. In a preferable application, the illuminance spectrum le(X) of the environmental light can be 
determined in advance as a prerequisite for color reproduction. In this case, unknown variables in Equation 
4 are only the angle of reflection e and the angle of deviation p, which will be determined by setting three- 
dimensional relative positions of the observer, the printed image, and the light source. Substitution of the 
angle of reflection e and the angle of deviation p thus determined into Equation 4 will give the illuminance 
spectrum 1(6, p.X) of the reflected light. 

With the illuminance spectrum 1(0, p,X) thus determined, tristimulus values X, Y, and Z in the CIE-XYZ 
colorimetric system are given by: 



10 



15 



rlSO 

X = kj Ue.pAV dk . • • ( 5a) 



^•780 

Y=k| I(e,pA)* y(MdX ...(5b) 
J380 



25 



Jr780 
I(e,p,X)- z(X)dX ...(5c) 
380 



30 



where x(X)., y(X), and z(X) denote color matching functions. 

Conversion of the tristimulus values X, Y, and Z obtained by Equations 5a through 5c, or their 
chromaticity coordinates (x,y,2), into color data of a colorimetric system, for example, the RGB colorimetric 
35 system, of an output device will precisely reproduce the colors of the printed image by the output device. 

B. Second Method of Color Computation 

The second method of color computation executes the integration of the first method in advance and 
40 utilizes the results of the integrations in computing colors. 

Substitution of Equation 4 into Equations 5a through 5c gives: 



45 



so 



55 



Bf4SDOCID: <EP_0669754A2J_> 



75 



20 



25 



30 



35 



40 



EP 0 669 754 A2 



Jr780 flBO 
I Ss(X)- <|)(X.)- x(X) dX + cos e I Sb(X)- x(X) dX 

380 hno 



J rim 
I Ic(X)- x(X) dX > ..•(6a) 
380 



i*780 

I Ss(X)- y(X) dX + cos 8 I Sb{X)- <t)(X)- y(X) dX 

380 /380 



Jrf-780 
I Ic(X)- y(X.) dX } ...(Sb) 
380 



Jr780 i-780 
1 Ss(X)- <|)(X)- z(X) dX + cos e I Sb(X)- <1)(X)- z(X) dX 
380 7380 



Jr780 
380 



780 

Ie(X)- Z(X) dX } ...(6c) 



The terms of the right hand side of each of Equations 6a through 6c can be rewritten with tristimulus 
values (Xb,Yb,Zb) of the internal reflection light, tristinnulus values (Xs,Ys,Zs) of the specular reflection light. 
45 and tristimulus values (Xe,Ye,Ze) of the environmental light as follows: 



50 
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r780 

(Xs,Ys,Zs) = k| Ss(X)- 4>(X)- (x(X),y(X).Z(X)) dk • . . (7a) 

/380 



..Ys.Z,) = k 



( 

5,Yi>,Zb) = k| 



r780 

(Xb,Yi>,Zb) = k| s^ixy {>r(X),y(X),z(>.)) dx ... (7b) 

'3 80 



i»7KU 

(Xe, Ye,Ze) = k [ Ie(X)- (x(X) ,y(X) .Z(?l)) dX . . . ( 7C ) 

/38O 



Using Equations 7a through 7c. the tristimulus values X, Y, and Z are given by: 



X = k{COS"p • Xs + COS0 • Xb + Xe} 
Y = k{cos"p • Ys + COS0 • Yb + Ye} 



(8a) 

Y = k{cos"p • Ys + COS0 • Yb + Ye} (8b) 

25 Z = k{COS"p*Zs + COS0»Zb + Ze} (8c) 

The sunn of the tristimulus values X + Y + Z is expressed as: 

X + Y + Z = k{cos"p(Xs + Ys + Zs) + cos0(Xb + Yb + Zb) + (X^ + Ye + Z^)} (9) 

30 

Since Equations 7a through 7c do not include the angle of reflection 6 and the angle of deviation p, the 
three sets of tristimulus values (Xs,Ys,Zs), (Xb.Yb.Zb), and (Xe,Ye,Ze) can be determined in advance by 
setting the coefficients of reflection Ss(X) and Sb(X). the luminance spectrum <^>(\) of the light source, and the 
illuminance spectrum le(X) of the environmental light and by executing the integrations of Equations 7a 
35 through 7c. 

The three sets of tristimulus values (Xs,Ys,Zs), (Xb,Yb,Zb). and (Xe,Ye,Ze) thus obtained are used to 
determine the tristimulus values (X,Y,Z) of the reflected light by substituting the angle of reflection 0 and the 
angle of deviation p into Equations 8a through 8c. Since the second method executes the integrations prior 
to the color computation, it has such an advantage that only the relatively simple operations are required for 
40 determination of the tristimulus values (X,Y,Z) of the reflected light in the actual color reproduction process. 

In the second method of color computation, chromaticity coordinates (Xg.yg), (Xb.yt). and (Xe,ye) may be 
used in place of the tristimulus values (Xs.Ys.Zs). (Xb,Yb,Zb), and (Xe.Ye.Ze). 



C. Third Method of Color Computation 

Based on the substitution given by the following Equations 10a through 10c. the right hand side of 
Equation 9 is rewritten as Equation 1 1 given below: 

Xs + Y3 + Zs = Ts (10a) 

50 Xb + Xb + Zb = Tb (10b) 

Xe + Ye + 2e = Te (10C) 

X + Y + Z = k{cos"p • Ts + cose • Tb + Te} (11) 

Equations 8a through 8c and Equation 11 give the chromaticity coordinate x corresponding to the 
55 tristimulus values (X.Y.Z) as: 
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X = X k{cos" p Xs + cos e -Xh X e} 

Y -I- V -I. '7 " ' • - . ( 12 ) 

A + Y + i^^cOS^ p + COS 0 Tb + Te} 

Equation 12 can be rewritten as Equation 13 by using the general relation specified as Equation 14 and 
by defining Cs, Cb, and Ce by Equations 15a through 15c: 



^ _ k{cos" p -Xs- Ts + cos e -Xb- Th + Xg- Tp} 
k{cos" p Ts + cos e -Tb + Tc} 

= Cs- Xs +Cb- Xh +Cc- Xe . . • ( 13 ) 



Xs = Xs"[-Kc-[Ts (14a) 
Xb = Xb-Tb (14b) 

Xe = Xe-Te (14c) 



Q ^ cos" p Ts 



cos" p Ts + cos 6 Tb + Te 



(15a) 



^ cos 0 • Tb 

^b- . . . (15b) 

cos" p Ts + cos e Tb + Te 



Ce = • ..(15C) 

40 cos" p Ts + cos 0 Tb + Tc 

The sinnilar equations to Equation 13 will be given for the other chromaticity coordinates y and z. The 
chromaticity coordinates (x,y) of the reflected light are thus given by: 



(x,y) = Cs' (Xs.ys) + Cb* (Xb.yb) + Ce* (Xe,ye) (16) 

Since x + y + z = 1 holds fronn the definitions of the chromaticity coordinates, the chronnaticity coordinate 
z is also obtainable from Equation 16. 
50 The meaning of Equation 16 will be easily understandable under the ideal observation condition without 

any environmental light. The third term of the right hand side of Equation 16 is equal to zero under the ideal 
observation condition with no environmental light. In this case, Equation 16 and Equations 15a through 15c 
can be rewritten respectively as Equation 17 and Equations 18a and 18b given below: 

55 (x,y) = Cs- (Xs.ys) + Cb* (Xb.yb) (1 7) 



10 
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10 



COS" p Ts ...{18a) 
cos" p -Ts + cos 6 Tb 



cos e - Tb ...(18b) 

cos" p -Ts + cos 6 Tb 



Equation 17 gives the chromaticity coordinates (x,y) of the reflected light in observation of the reflected fight 
with the vision system shown in Fig. 2 under the ideal condition with no environmental light. 

75 As clearly understandable from Equations 7a through 7c. (xs.y^) and (Xb,yb) of Equation 17 respectively 
represent the coordinates of the specular reflection light Ls and those of the internal reflection light Lb. Fig. 
3 is a chromaticity diagram of the chromaticity point Ps(Xs,ys) of the specular reflection light Ls and the 
same Pb(Xb.yb) of the internal reflection light Lb under the ideal condition with no environmental light. As 
illustrated in Fig. 3, the chromaticity point Pr(x,y) of the reflected light given by Equation 17 exists on a line 

20 drawn between the chromaticity point Ps(xs,ys) of the specular reflection light Ls and the chromaticity point 
Pb(Xb,yb) of the internal reflection light Lb. The chromaticity point Pr of the reflected light is the point of 
division between the chromaticity points Ps and Pb as described below. 

The coefficients of reflection Ss(X) and Sb(X) used in Equations 7a through 7c are determined by 
measuring the reflection factors with a light source having the predetermined luminance spectrum <>(x) 

25 according to the process described below under the ideal condition with no environmental light. The two 
sets of tristimulus values {Xs,Ys.Zs) and (Xb.Yb.Zb) are determined by substituting the coefficients of 
reflection Ss(X) and Sb{X) and the luminance spectrum <t>(\) into Equations 7a through 7c. The chromaticity 
coordinates (x^.y^) of the specular reflection light Ls and the chromaticity coordinates (Xb.yt) of the internal 
reflection light Lb included in the right hand side of Equation 17 can be determined according to the 

30 definitions of the tristimulus values {Xs,Ys.2s) and (Xb,Yb,Zb) given as Equations 14 and 10a through 10c. 

The coefficients Ts and Tb of the right hand side of Equations 18 are constant as shown by Equations 
10, and the coefficients Cs and Cb therefore depend upon only the angle of reflection 0 and the angle of 
deviation p. Since both the angles e and p range between 0 degree and 90 degrees, the coefficients Cs and 
Cb have always positive values. The sum of the coefficients Cs and Cb is equal to one under the condition 

35 with no environmental light as follows: 

Cs + Cb=1 09) 

The chromaticity point Pr(x,y) of the reflected light given by Equation 13 and the similar equation for y 
40 is the point of division between the chromaticity points Ps(Xs.ys) and Pb(Xb,yb) by the coefficients Cs and 
Cb. which are determined according to the angle of reflection e and the angle of deviation p. As described 
previously, the chromaticity point Ps of the specular reflection light Ls and the chromaticity point Pb of the 
internal reflection light Lb may be determined in advance through the measurement. Accordingly, the 
chromaticity coordinates (x.y) of the reflected light from a color print, which is disposed in the three- 
45 dimensional space, will be obtained by substituting the angle of reflection 0 and the angle of deviation p into 
Equations 17 and 18. 

In reproduction of a color print under practical observation conditions with environmental light, Equa- 
tions 16 and 15a through 15c are used instead of Equations 17 and 18. Upon condition that the coefficient 
Te of the environmental light is a constant, the coefficients Cs, Cb, and Ce are determined with the angle of 
50 reflection 6 and the angle of deviation p. The chromaticity coordinates (xe.ye) of the environmental light can 
be set to a specific value representing the environment. The chromaticity coordinates (x,y) of the reflected 
light will be thus determined according to the magnitudes of the angle of reflection d and the angle of 
deviation p. 

In one preferable application, plural sets of environmental light-based components Te, Xe, ye represent- 
55 ing typical observation environments such as print shops and offices are prepared in advance, and an 
appropriate one is selected among the plural sets in color reproduction. 

Since the second method and the third method of color computation do not require the integration 
operations, they can determine the chromaticity of the reflected light more readily and quickly than the first 
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method of color computation. 

Any of the three methods of color computation is applicable to determine the chromaticity of a color 
print disposed in the three-dimensional space from the angle of reflection e and the angle of deviation p of 
the observation direction from the reflection direction. Accordingly, the above methods can simulate the 
observation of the color print from various angles. When the angle of deviation p is substantially equal to 
zero degree, for example, the above methods can simulate the glossy condition of the color print where the 
component of the specular reflection light Ls becomes significantly large. 

As is generally known, a model called 'standard dichromatic reflection model' is applicable to reflection 
from such objects as plastics, pottery, and polyvinyl resins. The 'standard dichromatic reflection model' is 
constructed on the assumption that the specular reflection coefficient Ss in Equation 2 is constant 
irrespective of the wavelength X. The 'standard dichromatic reflection model' cannot, however, be applied to 
color prints, which require considering the dependence of the specular reflection coefficient Ss on the 
wavelength X. The first through the third methods of color computation described above take into account 
the wavelength dependence of the specular reflection coefficient Ss to compute the colors of a color print 
These methods accordingly allow the colors of the color print to be reproduced more precisely than the 
conventional method. 

D. Measurement of Reflection Properties of Printed Image 

For determination of the reflection coefficients Ss(X) and Sb(X) in Equations 12. print samples were 
prepared and their reflection properties were measured. The total of six samples were prepared by solidly 
pnntmg inks of yellow, cyan, and magenta on two different types of printing paper, that is, art paper and 
coated paper, respectively. 

Fig. 4 conceptually shows a measurement system for reflection properties. While a sample 2 placed on 
an optical sample table 1 was irradiated with light emitted from a light source 3, the illuminance or 
Illuminance spectrum of the reflected light was measured with a spectral radiance meter 4. The data 
obtained were analyzed by a personal computer 5. The light source 3 used here was a day-light flood light 
of standard illuminance D65. The measurement was conducted in a dark room to realize the ideal 
observation condition without any environmental light. 

The illuminance spectra of the reflected light were measured with the measurement system shown in 
Fig. 4 under a plurality of conditions that the angle of reflection 9 and the angle of deviation p are varied 
parametncally. Fig. 5 is a graph showing an example of the illuminance spectra obtained at three different 
values of the deviation angle p for the sample prepared by solidly printing cyan ink onto coated paper. As 
Illustrated in Fig. 5, the shape of the spectra varies with the angle of deviation p. The illuminance spectrum 
does not, on the other hand, change the shape in response to the change of the angle of reflection e but 
uniformly increases or decreases as the angle of reflection e. The illuminance spectra of Fig. 5 correspond 
to the Illuminance spectrum l(e.p,X) of the reflected light in Equation 3. 

Analysis of the data of illuminance spectrum 1(9, p,x) by Singular Value Decomposition (SVD) was then 
performed to determine whether the illuminance spectrum l(9,p,x) could be expressed by linear combination 
of variable-separated components as Equation 3. The results of the Singular Value Decomposition showed 
that the contribution rate up to the second singular value was extremely high as 0.9996 or greater in any 
one of the six samples measured. This proved that the illuminance spectrum l(fl,p,x) of the reflected light 
could be expressed in the form of Equation 3. 

The specular reflection coefficient Ss(X) and the internal reflection coefficient Sb(X) will be determined 
by solving simultaneous equations, which are obtained by substituting into Equation 3 plural sets of data of 
Illuminance spectrum l(e,p,x) collected under such condition that the angle of reflection e (or the angle of 
incidence e) is varied while the angle of deviation p is kept equal to zero degree. Alternatively, the specular 
reflection coefficient Ss(X) and the internal reflection coefficient Sb(X) can be determined by a method of 
least squares. In the latter method, Ss(X) and St,(\) having the least squares sum in Equation 3 are 
determined from the plural sets of the data of illuminance spectrum l(e,p,x) obtained under the condition 
that the angle of reflection $ is varied while the angle of deviation p is kept equal to zero degree. 

The exponent n of cosp will be determined by the least squares method from the plural sets of the data 
of Illuminance spectrum l(e,p,x) obtained under the condition that the angle of reflection 0 is kept constant 
while the angle of deviation p is varied. 

Figs. 6 and 7 are graphs showing the chromaticity data of the reflected light obtained for the six 
samples. Fig. 6 and Fig. 7 respectively show the results of coated paper samples and art paper samples. 

In the results of the coated paper samples shown in Fig. 6. the chromaticity data of the reflected light 
substantially follow a straight line for each ink of yellow, cyan, and magenta. The chromaticity point of the 
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light source is also shown In Fig. 6 for the reference. The terminal points of the chromaticity data of each 
ink correspond to the chromaticity point Ps of the specular reflection light and the chromaticity point Pb of 
the internal reflection light shown in Fig. 3. As is generally known, upon condition that the standard 
dichromatic reflection model having a constant specular reflection coefficient Ss is applicable, the chro- 

5 maticlty point Ps of the specular reflection light is identical with the chromaticity point of the light source. 
On the contrary, none of the straight lines, which the chromaticity data of the respective inks follow, goes 
through the chromaticity point of the light source In Fig. 6. This shows that the standard dichromatic 
reflection model Is not applicable to actual color prints on paper. 

The results of Fig. 6 clearly show that the data of chromaticity follow a straight line in each ink sample 

JO printed on coated paper and therefore that the third method described above is favorably applicable to 
reproduce the colors of a print precisely. This is also true for the art paper samples shown in Fig. 7. 

D. Color Reproduction Method 

75 Fig. 8 is a block diagram showing a computer system embodying the invention, which Is used for 

reproducing a color print disposed in the three-dimensional space. The computer system includes a CPU 
10 and a bus line 12. The bus line 12 is coupled to a ROM 14 and a RAM 16 and via input-output interfaces 
40 and 41 to a keyboard 30, a mouse 32, a digitizer 34, a color CRT 36, a magnetic disk 38. and a full color 
printer 42. The CPU 10 executes application programs stored in the ROM 14 or loaded to the RAM 16 to 

20 implement the processing described below. The computer system performs a rendering process to 
reproduce a color print disposed in the three-dimensional space. 

Fig. 9 shows a three-dimensional orientation of a color print 50, a light source 52, a screen 54 of the 
CRT, and an observation point 56 which is used in rendering by Ray Tracing. Ray Tracing Is a known 
rendering process In three-dimensional computer graphics, which determines color data (RGB data) of each 

25 pixel PX on the screen 54 by following a ray line RL running through each pixel PX on the screen 54 and 
the observation point 56 reversely from the observation point 56 to the light source 52, The print 50 Is 
previously divided Into a plurality of small polygons EL. Although each polygon EL is flat, the whole print 50 
may be modeled to have a three-dimensional curve. The angle of incidence 0 and the angle of deviation p 
are defined at the small polygon EL which exists on an Intersection of the ray line RL and the print 50 as 

30 shown in Fig. 9. 

Fig. 10 is a flowchart showing a procedure for computing color data (RGB data) of each pixel PX on the 
screen 54 by Ray Tracing. At step SI, one of the three methods of color computation is selected. 

When the process according to Equation 4 is selected, the program goes to step S2 at which the 
reflection coefficients Ss(X) and Sb(X) are read out from the RAM 16 according to the angle of reflection e 
35 and the angle of deviation p, or interpolated by the angle of reflection e and the angle of deviation p, and 
the illuminance spectrum l(0,p,X) of the reflected light is determined. Data used in Equation 4, that is, the 
reflection coefficients Ss(X) and Sb(X), the exponent n, the luminance spectrum <^(X) of the light source, and 
the illuminance spectrum le(X) of the environmental light, are loaded in advance from the magnetic disk 38 
to the RAM 16. 

40 At step S3, tristimulus values X, Y, and Z are determined by executing the integration of the illuminance 
spectrum l(0,p,X) by Equation 5. The program then proceeds to step SB at which the tristimulus values X, Y, 
and Z are converted to RGB data and allocated to a pixel on the screen 54 of the CRT. In a concrete 
operation, RGB data are stored at each pixel position in a frame memory. The conversion from the XYZ 
colorimetric system to the RGB colorimetric system depends upon the properties of the output device. 

45 Application of a conversion equation suitable for the output device allows the precise color reproduction in 
the specific color reproduction range of the output device. 

When the process according to Equations 8a-8c is selected, the program goes to step S4 at which the 
tristimulus values (X.Y,Z) of the reflected tight are determined by reading out the tristimulus values 
(Xs.Ys.Zs), (Xb.Yb.Zb), and (Xe.Ye,Ze) from the RAM 16 according to the angle of reflection 0 and the angle 

50 of deviation p, or by interpolating the tristimulus values (Xs.Ys.Zs). (Xb.Yb.Zb), and (Xe,Ye,Ze) by the angle 
of reflection 6 and the angle of deviation p. Data used in Equations 5a-5c, that is, the tristimulus values 
(Xs,Ys,Zs). (Xb.Yb.Zb). and (Xe,Ye,Ze), the exponent n, and the factor k, are loaded in advance from the 
magnetic disk 38 to the RAM 16. The program then proceeds to step SB at which the tristimulus values X, 
Y, and Z are converted to RGB data and allocated to a pixel on the screen 54 of the CRT. 

55 When the process according to Equations 15a-15c and 16 is selected, the program goes to step 85 at 
which the chromaticity coordinates (x.y) of the reflected light are determined by reading out the chromaticity 
coordinates (Xs.ys), (Xb.yb). and (Xe.ye) from the RAM 16 according to the angle of reflection e and the angle 
of deviation p. or by interpolating the chromaticity coordinates (Xs.ys). (Xb.yb). and (Xe,ye) by the angle of 
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reflection e and the angle of deviation p. Data used in Equations I5a-I5c and 16, that is, the chromaticity 
coordinates (Xg.ys), (Xb.Vb). and (Xe,ye), the exponent n, and the factors Ts, Tb. and Te, are loaded in 
advance from the magnetic disk 38 to the RAM 16. The program then proceeds to step S6 at which the 
chromaticity coordinates (x.y) thus obtained are converted to RGB data and allocated to a pixel on the 
5 screen 54 of the CRT. In the process according to Equations 15a- 15c and 16. the chromaticity coordinates 
(x,y) determine only the relative values of RGB data but no absolute values. Accordingly, the absolute levels 
of RGB data will be determined separately. 

As described above, the tristlmulus values X. Y, Z in the device-independent CIE-XYZ colorimetric 
system or their chromaticity coordinates (x.y) can be obtained by any of the three methods of color 
10 computation. By converting the tristimulus values X. Y. Z or the chromaticity coordinates (x.y) to another 
colorimetric system (RGB colorimetric system, YMCK colorimetric system, or the like) with a conversion 
equation suitable to the output device, such as a color CRT and a printer, the color of the color print is 
precisely reproduced with the output device within its specific color reproduction range. 

The above embodiment is only illustrative and not restrictive In any sense. There may be many 
15 modifications, alterations, and changes without departing from the scope or spirit of essential characteristics 
of the invention, for example, as follows: 

(1) Any other device-independent colorimetric systems such as CIE-L'aV colorimetric system and CIE- 
L*uV' colorimetric system are applicable in place of the CIE-XYZ colorimetric system used in the above 
embodiment. 

20 (2) The reflection coefficients Ss(X) and Sb(X) can be determined by the following method. Illuminance is 
first measured for each color Ink, that is, yellow, cyan, and magenta. The chromaticity coordinates (x.y) 
are then calculated from the measured illuminance and plotted in an x-y chromaticity diagram as shown 
in Fig. 6 or 7. The inner-most chromaticity point in the series of chromaticity coordinates (x,y) Is selected 
and the illuminance spectrum of the selected point is used to determine a spectral reflection coefficient, 

25 which is used as the specular reflection coefficient Ss(X). The outer-most chromaticity point in the series 
of chromaticity coordinates (x,y) is selected and the illuminance spectrum of the selected point is used 
to determine another spectral reflection coefficient, which is used as the internal reflection coefficient Sb- 
(X). These reflection coefficients Ss(X) and Sb(X) are determined for each color ink. accordingly. This 
method realizes simple determination of the reflection coefficients Ss(X) and Sb(X). 

30 (3) Although only one light source 52 is used in the color reproduction process of the embodiment, a 
plurality of light sources may be used for computation of color data of each pixel PX on the screen 54 to 
attain more natural color reproduction. In the actual room or natural environment, the condition that the 
print is irradiated with light emitted from only one light source is seldom found, but the print Is generally 
Irradiated with various light. In the color reproduction process with a plurality of light sources, color data 

35 of each pixel PX are obtained for each light source according to the routine shown in Fig. 10 and then 
accumulated. 

Although the present Invention has been described and illustrated in detail, it is clearly understood that 
the same is by way of Illustration and example only and is not to be taken by way of limitation, the spirit 
and scope of the present invention being limited only by the terms of the appended claims. 
40 The features disclosed in the foregoing description, in the claims and/or in the accompanying drawings 
may, both separately and in any combination thereof, be material for realising the invention In diverse forms 
thereof. 

Claims 

45 

1, A method of simulating a color of each reflection point in a color print observed by an observer with an 
output device, the color print being disposed in a three-dimensional space and irradiated by a light 
source, the light source emitting a light beam having a luminance spectrum <t>{\) where X denotes a 
wavelength of the light beam, the method comprising the steps of: 
50 determining an illuminance spectrum l(a,p,x) of reflected light reflected at the reflection point 

according to an Equation: 

l(0.p.X) = {Ss(X)* COS"p + Sb(X).COS0}<|>(X) + le(X) 

55 wherein Sb(X) and Ss(X) denote first and second reflection coefficients, respectively, p denotes an angle 

between a first direction of the reflected light and a second direction from the reflection point to the 
observer, n is a constant, d denotes an angle of reflection of the reflected light, and le(X) denotes an 
illuminance spectrum of environmental light observed by the observer; 
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integrating the illuminance spectrum l(9.p.X) over a wavelength range of visible light with respective 
three color matching functions of a device-independent colorimetric system, to thereby obtain 
tristimulus values in the device-independent colorimetric system representing the color of the reflection 

point; and 

5 producing an image of the color print with the output device as a function of the tristimulus values. 

2. A method in accordance with claim 1 . wherein 

the first reflection coefficient Ss(X) regarding specular reflection light depends upon the wavelength 

X. 

10 

3. A method in accordance with claim 1, further comprising the step of: 

converting the tristimulus values to color data in a specific colorimetric system of the output device. 

4. A method in accordance with claim 1 , wherein 

15 the illuminance spectrum le(X) of environmental light is equal to zero. 

5. A method of simulating a color of each reflection point in a color print observed by an observer with an 
output device, the color print being disposed in a three-dimensional space and irradiated by a light 
source, the light source emitting a light beam having a luminance spectrum <^>(X) where X denotes a 

20 wavelength of the light beam, the method comprising the steps of: 

determining target tristimulus values X, Y, and Z of a device-independent colorimetric system for 
reflected light reflected at the reflection point according to Equations: 

X = k{cos"p • Xs + cose • Xb + Xe} 
25 Y = k{cos"p*Ys + cos0»Yb+ Ye} 

Z = k{COS"p • Zs + COS0 .Zb = Ze} 

where (Xs,Ys,Zs), (Xb,Yb,Zb). and (Xe,Ye.Ze) respectively denote first, second, and third sets of tristimulus 
values, k is a first constant, p denotes an angle between a first direction of the reflected light and a 
30 second direction from the reflection point to the observer, n is a second constant, and d denotes an 
angle of reflection of the reflected light; and 

producing an image of the color print with the output device as a function of the target tristimulus 
values X, Y. and Z. 



35 6. A method in accordance with claim 5, further comprising the step of: 

converting the target tristimulus values X, Y, and Z to color data of a specific colorimetric system of 
the output device. 

7. A method in accordance with claim 5, wherein 

40 the third set of tristimulus values (Xe,Ye,Ze) regarding environmental light are equal to zero. 

8. A method of simulating a color of each reflection point in a color print observed by an observer with an 
output device, the color print being disposed In a three-dimensional space and irradiated by a light 
source, the light source emitting a light beam having a luminance spectrum ^(X) where X denotes a 

45 wavelength of the light beam, the method comprising the steps of: 

determining target chromaticity coordinates (x.y) of a device-independent colorimetric system for 
reflected light reflected at the reflection point according to Equations: 



(x.y) = Cs- (Xs.ys) + Cb* (Xb.Yb) + Ce*(Xe.ye) 

50 
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^ . COS" p Ts 





COS" p Ts + COS e Tb + T, 



e 
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15 



COS 9 • Tb 

cos" p Ts + cos e -Tb + Te 



Ce = 

COS" p Ts + COS e Tb + Te 



where (Xs.ys), (Xb.yb), and (xe.ye) denote first, second, and third chromaticity coordinates. Tg, Tb, and Tg 
are constants, p denotes an angle between a first direction of the reflected light and a second direction 
20 from the reflection point to the observer, n is a constant, and e denotes an angle of reflection of the 
reflected light: and 

producing an innage of the color print with the output device as a function of the target chromaticity 
coordinates (x.y), 

25 9. A method in accordance with claim 8, further comprising the step of: 

converting the target chromaticity coordinates (x,y) to color data of a specific colorimetric system of 
the output device. 

10. A method in accordance with claim 8, wherein 

the third set of chromaticity coordinates (Xe,ye) and the constant Te regarding environmental light 
are equal to zero. 

11. A method in accordance with claim 1, wherein 
the device-independent colorimetric system is CIE-XYZ colorimetric system. 



35 



12. 



40 



An apparatus for simulating a color of each reflection point in a color print observed by an observer, the 
color print being disposed in a three-dimensional space and irradiated by a light source, the light 
source emitting a light beam having a luminance spectrum 0(X) where X denotes a wavelength of the 
light beam, the apparatus comprising: 

means for determining an illuminance spectrum l(0,p,X) of reflected light reflected at the reflection 
point according to an Equation: 

\{d,p,\) = {Ss(X)*COS"p + Sb(X).COS0}<^>(X) + le(X) 

45 wherein Sb(X) and Ss(X) denote first and second reflection coefficients, respectively, p denotes an angle 

between a first direction of the reflected light and a second direction from the reflection point to the 
observer, n is a constant, $ denotes an angle of reflection, and le(X) denotes an illuminance spectrum of 
environmental light observed by the observer; 

means for integrating the illuminance spectrum 1(0, p,X) of the reflected light over a wavelength 

60 range of visible light with three color matching functions, respectively, of a device-independent 

colorimetric system, to thereby obtain tristimulus values in the device-independent colorimetric system; 
and 

output means for producing an image of the color print as a function of the tristimulus values. 

55 13. An apparatus in accordance with claim 12, wherein 

the first reflection coefficient Ss(X) regarding specular reflection light depends upon the wavelength 
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14. An apparatus in accordance with clainr^ 12, further connprising: 

means for converting the tristimulus values to color data of a specific colorimetric system of said 
output means. 

5 15. An apparatus in accordance with claim 12, wherein 

the illuminance spectrum le(X) of environmental light is equal to zero. 

16. An apparatus for simulating a color of each reflection point in a color print observed by an observer, the 
color print being disposed in a three-dimensional space and irradiated by a light source, the light 
;o source emitting a light beam having a luminance spectrum <^(X) where X denotes a wavelength of the 
light beam, the apparatus comprising: 

means for determining target tristimulus values X, Y, and Z of a device-independent colorimetric 
system for reflected light reflected at the reflection point according to Equations: 



where (Xs,Ys,Zs), (Xb.Yb.Zb). and (Xe,Ye,Ze) respectively denote first, second, and third sets of 
20 tristimulus values, k is a first constant, p denotes an angle between a first direction of the reflected light 

and a second direction from the reflection point to the observer, n is a second constant, and e denotes 
an angle of reflection of the reflected light; and 

output means for producing an image of the color print as a function of the target tristimulus values 
X, Y, and Z. 

25 

17. An apparatus in accordance with claim 16, further comprising: 

means for converting the target tristimulus values X, Y, and Z to color data of a specific 
colorimetric system of said output means. 

30 18. An apparatus in accordance with claim 16, wherein 

the third set of tristimulus values (Xe,Ye,Ze) regarding environmental light are equal to zero. 

19. An apparatus for simulating a color of each reflection point in a color print observed by an observer, the 
color print being disposed in a three-dimensional space and irradiated by a light source, the light 
35 source emitting a light beam having a luminance spectrum ^(X) where X denotes a wavelength of the 
light beam, the apparatus comprising: 

means for determining target chromaticity coordinates (x,y) of a device-independent colorimetric 
system for reflected light according to Equations: 



75 



X = k{COS"p • Xs + COS0 • Xb + Xe} 

Y = k{cos"p • Ys + cose • Yb + Ye} 
Z = k{cos"p • Zs + cose • Zb + Ze} 



40 



(x,y) = Cs« (Xs.ys) + Cb- {Xb,yb) + C^* (Xe.ye] 



Cs = 



COS" p Ts 



45 



cos" p Ts + COS 6 Tb + Tc 



50 



Cb = 



COS B • Tb 



COS" p -Ts -t- COS e Tb + Tc 
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COS" p Ts + COS 6 Tb + Te 



17 



BNSOOCIO: <EP_0669754A2J^ 



EP 0 669 754 A2 



• 



where (Xs,ys), (Xb.Vb). and {Xe,ye) denote first, second, and third chromatlcity coordinates. T^. Tb. and Te 
are constants, p denotes an angle between a first direction of the reflected light and a second direction 
from the reflection point to the observer, n is a constant, and e denotes an angle of reflection of the 
reflected light; and 

s output means for producing an image of the color print as a function of the target chromaticity 

coordinates (x,y). 

20. An apparatus in accordance with claim 19, further comprising: 

converting the target chromaticity coordinates (x,y) in the device-independent colorimetric system 
10 to color data of a specific colorimetric system of said output means. 

21. An apparatus in accordance with claim 19, wherein 

the third set of chromaticity coordinates (Xe,ye) and the constant Tg regarding environmental light 
are equal to zero. 

15 

22. An apparatus in accordance with claim 12, wherein 

the device-independent colorimetric system is CIE-XYZ colorimetric system. 
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@ Method and apparatus for simulating color print 



@ There is provided a nnodel where the illuminance 
of reflected light fronn the print is expressed with a 
linear combination of the illuminance of specular 
reflection light and that of internal reflection light. On 
the basis of the model, three methods are applicable 
^ to compute the colors of a print disposed in the 
^ three-dimensional space. In the first method, a 
specular reflection coefficient and an internal reflec- 
tion coefficient, which are depending upon the 
1^ wavelength, are interpolated by an angle of reflection 
$ and an angle of deviation p, and the illuminance 
2 spectrum of the reflected light is subsequently deter- 
^ mined (step 82). Tristimulus values X, Y, and Z are 
then determined by integration of the illuminance 
spectrum according to the color matching functions 
^ (step S3). In the second method, the tristimulus 
^ values are determined for the specular reflection 
light, the internal reflection light, and the environ- 
mental light, respectively, and then interpolated by 



the reflection angle 0 and the deviation angle p (step 
84). In the third method, the chromaticity coordi- 
nates are determined respectively for the specular 
reflection light, the internal reflection light, and the 
environmental light, and then interpolated by the 
reflection angle 6 and the deviation angle p (step 
S5). 
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